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1 Introduction 

The study of strongly interacting matter at high temperature and density 
had attracted the interest of people all around the globe over the past couple 
of years. Now with the availability of the new experimental facilities such as 
the ep collider of HERA at DESY, Fermi-Lab Tevatron ( pp, p A collision ), 
BNL Relativistic Heavy Ion Collider( RHIC) and LHC of CERN; there is a 
renewed interest in the field-from experimental as well as theoretical point of 
view to observe as well as theoreticaly model the production of Quark Gluon 
Plasma -the elusive state of strongly interacting matter- that is supposed 
to exist in a chiral symmetry restored and deconfined phase at sufficiently 
high temperature or density or both together. So it is expected that with the 
availability of the new experimental data, it will be possible to pin down 
the uncertainties that existed in the previous estimates and hence provide a 
better model of the physics involved. The interest in these studies stem from 
both physics as well as astro-physics point of view, e.g to understand the 
structure of the interior of a neutron star, millisecond pulsar or for that matter 
in the early universe cosmology, when the primordial plasma was believed 
to exist in the form of quark gluon soup all around the universe. A study 
of this phase of matter in laboratory conditions is expected to shed some 
light on these issues and hence to the understanding of the universe around 
us. But to observe or detect the QCD plasma produced in the laboratory 
one needs to have a proper estimate of the energy density per unit voloume, 
initial temperature , entropy density and the number density of the produced 
partons Q. More over the further space time evolution of the plasma as well 
as the experimental signatures also depends on these quantities. Therefore it 
is crucial to have a proper understanding of the Pre-equilibrium production 
phase of the QCD plasma. The essential goal of these studies is to understand 
how in a nucleus nucleus collision, the coherent distribution of partons in each 
nucleus evolve to a highly incoherent state, and finally through self interaction 
generates a collective QCD excitation -that with the passage of time evolves 
and re-hadronizes before getting captured in the detectors. Generally, the pre- 
equilibrium phase is supposed to start from the time of the nuclear contact 
and last till the system has thermalised. As has already been mentioned, all 
the planned experiments are expected to cover different energy scales ranging 
from 200 GeV to lOOTev, so as one moves from the SPS to LHC energy scale 
different paradigms of QCD opens up, from non-perturbative to perturbative! 
Since the time taken by the nuclei to pass through each other in LHC scale ~ 
0.0005 fm/c and the time of thermalisation( for gluons) ia around 0.25 fm/c, 
the system spends an appreciable amount of time in the pre-equilibrium phase. 
Therefore, we in our "qualitative" description of the pre-equilibrium phase of 
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QGP — will start from the low energy scale where the non-perturbative QCD 
mechanisms are dominant and move towards the high energy scale where the 
pQCD based processes dominates the scenario. 



1.1 Non-Perturbative Domain: The Flux Tube Model 

The idea of creating QGP throgh flux tube decay owes its existence to Low 
and Nussinov [0]. According to this model in a relativistic heavy ion colli- 
sion, as the target and the projectile nucleui passes through each other, due 
to multiple exchange of soft gluons they get color charged-creating a strong 
chromo-electric field between them. This field pulls out q, q pairs from the vac- 
uum by tunneling. A simple kinematic consideration shows, to create a pair of 
particles from vacuum the amount of energy required is~ 2m on the other hand 
as these charged particles move apart over a characteristic length scale, say 
Compton wave length, the energy gained by the system is ~ 2E/m. Now one 
can deduce from a simple energy conservation argument that critical strength 
field strength required to produce pairs from vacuum is E c ~ 2m 2 ||.In 
technical language it is said-if the strength of this electric field gE exceeds a 
critical value gE c = 2m 2 , the vacuum becomes unstable against production of 
q q pairs of mass m each. Since the system looks like a parallel plate capacitor 
,this model is also known as capacitor plate model. In this model the electric 
field is a c number field and assumed to be constant all through out. The rate 
of pair production per unit volume and per unit time is, 
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On the other hand because of self interacting nature of the non-Abelian gauge 
fields-the color-electric field set up between the receding nuclei can oscillate 
in time. In fact in this model, the surface charge density on the nuclear 
plates, in the transverse directions are imagined to be homogeneous, with 
plates extended up to infinity. As the plates move apart, one can find out the 
field configuration in the central Baryon free region, by solving the source- 
less Yang-Mills (Y.M) equations. In an idealised situation one can study the 
system in (1 + 1) dimension with an axial gauge choice to get, 
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Here since u = J (3\ 2 + 02 2 + Ps 2 , the frequency of oscillation is seen to be 
dependent on amplitude. This is the motivation to re-investigate Schwinger 
mechanism in an oscillating field. Since finding out W p using these exact 
solution, is a formidable task an analysis was carried out in [fjj, to see effect 
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of a sinusoidally time varying field on pair production. The result showed that, 
the production probability-due to the presence of the oscillating field goes like 

2m 

(3) 

Since u)q is the energy associated with each quanta of the oscillating gluonic 
field, so m divided by wo is the number of gluon quanta required to produce 
a pair, hence it is called multi-gluon ionization process. This result is valid 
even when E <C E c and dominates over that of the tunneling rate. A compar- 
ison with the perturbative approach also established that for constituent 
quark mass the multi Gluon ionisation process is dominant over the other two. 
The effect of time variation of the electric field was also recently considered 
in ||, for an exponentially decaying Electric filed ( Decay due to depletion of 
the electric field caused by particle production) and the conclusion was, time 
dependent field increases the production rate. We will discuss this in an ap- 
propriate section. Before closing the discussion on nonperturbative production 
mechanisms its worth mentioning that using similar techniques, one can es- 
timate the decay rate for Gluons too and for an constant external electric 
field, the rate was found to be H similar to the quark case with mass set to 
equal to zero and a different color factor sitting in front. Few other interest- 
ing extensions of this model have been performed inj|] where the effects of 
finite size correction and the effect of a confinning potential have been studied 
separately. 
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2 Hard and semi Hard processes 

Having outlined the non-perturbative production mechanisms, we will discuss 
here the perturbative mechanisms, responsible for the production of quarks 
and gluons at a higher energy scale. The idea here is, right after collision 
the coherent distribution of the partons over the nuclei gets disrupted, as 
the valance partons of each individual nuclei starts colliding with that of the 
other-producing an incoherent bunch of partonic shower that-decays either 
forming space like or time like cascades, producing a mini-jetty enviornment . 
The important input in this theory comes from the nucleonic structure func- 
tions measured or predicted for different values of the fractional longitudinal 
momentum that the colliding parton in question is carrying, as known in the 
literature as Bjorken x and the energy scale of collision Q 2 . The quantity of 
interest in this model is the jet cross section defined in terms of the double 
differential cross section for parton parton collisions ||, as, 

a ptpd y\d 1/2 at 

Here the subscripts a and b refer to the different species of partons. The 
Bjorken variable x\ and X2 denotes the longitudinal momentum fraction car- 
ried by the partons of respective nucleon with f a and fb as their structure 
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functions. Also — is the double-differential cross section associated with 
partons of type a and b. The energy density is estimated from eqn.(§) by 
using the relation; 

,_ dE T A _TAA(b) 

With TAA(b) as the impact parameter and aj e t(Ej?) energy moment of the 
cross section. The parton parton subprocess cross section used here has a 
pathological divergence, that goes as ~ Pt~ 4 - In order to get rid of this sit- 
uation a "saturation criteria" (|Tcj| is usually invoked, that gives po — 2GeV 
for the LHC and po ~ lGeV forRHIC for Pb on Pb collision at zero impact 
parameter. This model predicts that the number density of gluons will out 
number the quarks. The mean transverse energy of gluons found to be ~3 GeV 
with the initial temperature ~ 1.1 GeV. 

2.1 McLerran Venugopalan Model 

Other than these two complementary processes of QGP production at RHIC 
and LHC energy scales, there also exists a third scenario that recently was 
pointed out by L.McLarren and Raju VenuGopalan in a series of papers pd| ], 
where it was argued that there is also another source of mid rapidity gluons 
in the central rapidity region of large "A" nucleus nucleus collision. Their 
argument was, as the surface density of mean charge squared fluctuation i.e fi 2 
goes as, fi 2 = 1.1 A3 fm~ 2 then for a sufficiently large A nucleus, Aq^ <C [i 2 

such that a s (fi 2 ) <C 1 so that the method of weak coupling expansion would 
still be valid. In this scenario the soft gluon fields produced by the color charges 
is classical and can be obtained by solving the classical Y.M equations and 
from there they compute the distribution function given by, 



L(fc T ,A)~Jr> UE (7) 



where 



where A is a cutoff scale associated with dynamical screening effects in dense 
partonic medium. In the appropriate kinematical domain the predictions of 




this model has been matched with that of pQCD results [12|. There has also 
been some effort to calculate the same using the coherent state description 
instead of the classical radiation formula, and they were shown to agree with 
each other [13] 
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3 Space Time Evolution 



Once there are sufficient number of partons in the system, they will interact 
with each other and undergo a space time evolution. This space time evolu- 
tion of the QCD plas ma, should be described by a gauge covariant transport 
equations. Following ( pi} ) one can define the gluon and quark distribution 
functions as. 
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And use the equations of motion to arrive at the gauge covariant kinetic equa- 
tions of Elze, Gyulassy and Vasak.From there one can arrive at the classical 
non abelian transport equation, after taking an ensemble average of the oper- 
ator valued equation and then setting terms of 0(h = 0):, 



pfDJ? (x,p) + g/2p»d/ [F^, W (x,p)] + = 



(10) 



Here [, ] + means anticommutator. The distribution function W(x,p) , is a 
hermitian matrix in color space and for the SU(2) case can be written in 
terms of a color singlet and triplet components as, 
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Here (G) = TrW (x,p) and (G a ) = Tr [X a W {x,p)\. Using Eq. (|TT[) one can 
write Eq. ( [TQ| ) in terms of a set of coupled partial differential equations. The 
coupling is oetween the singlet and triplet distribution function for quarks. 
They are as follows, 



p^(G)+gp fM d/[F^ a (G a )} = 
p^{G a ) + ge abc p»A\(G c ) + g Pfi d/ [2F^ a {G a )\ = 



(12) 



To take care of the the effect of collisions along with non perturbative produc- 
tion mechanism one needs to put a source and a collision term on the right 
hand side of this equation. 



p^D^W (x,p) + g/2p»d/ [F^, W (x,p)] + = S ource + C olli 



(13) 



Lots of studies have been performed in the past ( see for instance references 
[fl5f to (lTj ) to investigate the space time evolution of the plasma including the 
effects of source and the collision term or just the source. The popular form of 
the source used in the colorflux tube model is given as, 
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And the collision term is usually written [17] with the help of the equilibrium 
distribution function Wq(p) and the mean collision time tq as:, 
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In all these investigations it was observed that due to production of the par- 
ticles the external field shows an exponentially damped oscillation in proper 
time. Prom the numerical solution of those equations one could see the fea- 
ture of conversion of fields to particles and vice versa. To take care of this 
time varying external field, an artificial time dependence is introduced in the 
(Schwinger inspired) source term. Since in an exponentially decaying external 
fiel d, t he tunneling picture of pair production is not appropriate, the authors 
of p8| taking this into account derived a modified source term, using pertur- 
bation theory ( with a decay constant of 0.1 fm). On the other hand if the 
produced particle current is strong, it would react back on the system. To 
take care of the effect of quantum back reaction by the particle current a se- 
ries of investigations were performed by F. cooper, E. Mottola et al. (jn§). And 
there it was pointed out, that the incorporation of the Pauli blocking in the 
Schwinger inspired source term almost reproduces the exact field theoretical 
treatment. The pauli blocked source term they choose was of the form 
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The result of their |19|] exact calculation is shown in Figure 1. and 2. where 
the solid line represents the result of the exact field theoretic calculation and 
the dotted line shows the result of incorporating a Pauli-Blocked Source term 
(Eq. ( |l6| ) )to the Boltzman Vlasov equation. The half life of the field can 
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Figure 1: Current vs time 



Figure 2: Electric Field vs time 
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be calculated from Fig. 2, and the value ( for a constituent quark mass) 
comes close to that estimated in ( f2~I|| ), i.e ~ 4 to 5 fm. On the other hand 
for current quark mass, the same number comes out to be 240fm! After 
comparing these numbers it seems that the use of constituent quark mass 
is more reasonable for this model than the current quark mass. The current in 
fig 1 is seen to execute an oscillation in proper-time. The origin of this strong 
oscillation is the perhaps due to the fact that, in presence of an external field 
the negative and positively charged particles tends to move away in opposite 
directions hence creating a field with an opposite polarity. As the strength of 
this field, overcomes that of the external field these charges moves in reverse 
direction and this prosess keeps on continuing generating an oscillation in 
time. The peak initial temperature attained, for constituent quark masses, in 
this model seems to be around 800 MeV, i.e well above the critical energy 
density required to produce the plasma. 



3.1 Tunneling At Finite Temperature 

In the previous section we have discussed how the field decays with time. We 
also had noted from the study of Cooper et. al [19], that by the time the field 

has decayed to ^th fraction of its original value, the system has attained some 
temperature. Though this temperature is an oscillating function of proper 
time, still the time scale of oscillation is much larger that the time scale of 
production. In ref.([^lj),it was argued from the respective time scale analysis 
of the relevant processes e.g production, depletion and thermalisation that, by 
the time the system thermalises, some residual mean field will still be present 
there in the system. The existence of such mean field for a fermionic as well 
as a bosonic system has also been noted in the studies of Blaizot and Iancu 
p0 |. To see the effect of temperature, on Schwinger's tunneling mechanism, an 
investigation was carried out in (2]J , assuming the external filed to be constant 
and homogeneous every where; the effect of finite size corrections was also 
taken into account in [22 [.For cyllindricaly symmetric system of radius R, the 
tunneling rate at a distance p from the axis was found out to be:, 
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The Fig. 2a shows, how the rate varies with temperature and Fig. 2b shows 
how it varies as one moves radially outward. 




3.2 Color Oscillation 

In the limit of very low number density and very weak external field, one 
can study the collective oscillation of the plasma, by taking the momentum 



moments of equations (12) and from there generating a set of color hydro- 
dynamic equations. Before we go for generating the hydrodynamic equations, 
we normalize the 4-velocities for quarks in the following way 

rrarj n - + > rra JV^ga (x,p)d A p J p»G (x,p) d 4 p 

UU » = 1 Wlth mU = JG.(x, P )^ = IG(*,p)#p (18> 

After taking the zeroth moment of eqn.([l^) we arrive at 

d tl [G"}=0 and d»[G a *\+ge ab(i A ll h G* = Q (19) 
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Following Stewart [24] if one decomposes the 4- vectors G^ = mnU^ and G a 



mn a U^ one arrives at the mass and color continuity equations, 

d^[nU»}=0 and d„ [U»n a ] - ge abc A*n c U» = (20) 

Further by defining the ratio of the two charge densities as Q a = 2* and using 
eq.(^) one can obtain the color evolution equation , namely 

[Q a ] - ge^A^QJJ* = (21) 



Similarly from the first momentum moment of equation (12) one arrives at 



U*dU» = i-Fa^ria = ^-Fr^ a (22) 
m m 

Here repeated indices are summed up and J^ a = n\J tl Q a . These single 
flavour equations are very similar to the hydrodynamic equations of Kajantie 



and Montonen [25] and they can easily be generalised to the multiflavor case 
for quarks. 

3.3 Gluon Hydro-dynamical Equations 

To derive the gauge covariant transport equations for Gluon one follows, the 
same procedure as before-starting from eqn( ^). Upon using the equation 
of motion, and setting terms 0(h) to zero one arrives at the classical kinetic 
equation of the gluons,i.e, 



p^.D^G^ + |p ct c£ [ T UT , ] + = g ( T^G° V - G^TZ ) (23) 



Here [, ]+ means anti-commutator and = — ig[A^,} where, A a ^ = 

~~ ^abc^w' ~^uv = — */abc^A»^ an( ^ ^abc * s ^ ne antisymmetric structure constant 
for SU(2) and g is the coupling constant. In order to simplify things a bit 
further, one can assume the quantity G^ u to be symmetric w.r.t the Lorenz 
indices, which is technically known as spin equilibration ansatz, to write 

Gnv(x,p) =p ft p v G(x,p), (24) 

where G(x,p) is a Lorentz scalar, but is a 3 x 3 matrix in color space. One 
can chec k th at with the ansatz of spin equilibration, the right hand side of 
equation[p3] vanishes and the Gluon kinetic equation in component (color) 
notation takes the form : 

jfduG™ + g^.A% [f cma G an - G ma f can ] + i\p a dl [f^G an + / ean G ma ] 

xF CT e T = (25) 



2 By multiplying equation (|22|) by U" one can show U^U^ = constant;Norm&lis'mg the 
U's appropriately Eqn.([l8|) is satisfied. 
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In the above equation all repeated indices are summed over. From equation 
( B4j) we define the diagonal, antisymmetric and symmetric components as 
follows : 
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We want to emphasize here that the most general the hydrodynamic equations 
for Gluons would not be the same as the quarks, it is an artifact of the spin 
equilibration ansatz that they seem to come out to be the same here. In order 
to make some more progress now we will further assume that all the symmetric 
combinations of the distribution function are zero. With this assumption one 
can show that the equations (|25|) reduce to : 
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The set of equations for studying the hydrodynamic evolution of the plasma 
can be obtained as in the quark case after taking moments of the one particle 
distribution function as before with the only exception that for the massless 
gluons one will have:, 



U"U„ 



0. 



(29) 



Instead of going into the details of the derivation, we shall write down here 
final set of Gluon hydrodynamic equations:, The continuity equations are, 



d„ [nW] = [Q a ] - geabcA^QcU^ = 

And the force balance equation is 

(J 



E 



(30) 
(31) 



For stud ying the collective behavior of the system one needs to solve Eq.( 
and Eq.( |3l|) along with the Yang Mills equation 
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self consistently. The basic idea here is that due to the self interaction the 
high momentum modes generate a low momentum long wavelength mean field 
which in turn acts back as a source term for a mean Yang-Mills field equa- 
tion. In order to study the collective properties one needs to solve these equa- 
tions self consistently, to describes the proper evolution of the system. We 
have tried to solve these equations numerically, in 2+1 dimension. Initially 
susch a study was performed in [^6| in the stationary frame of the plasma 
in (1+1) dimension, and it showed^depending upon the initial conditions — a 
rich variety of regular as well as chaotic solutions. The motivation of going 
to one extra dimension was to understand how the system behaves under a 
small perturbation in one direction with an objective to understand the role 
of chaotic dynamics in equilibrating the system. 



Transverse Oscillation 

Longitudinal Oscillation 




,, „, nl _ Time In Units Or Plasma Frequency 

Time in Units Of Plasma Frequncy * ' 



Figure 3: Oscillation Profile:Left panel V x vs time. Right panel V z vs time 



4 Study of Collective Oscillation of the Plasma 

The figures five and six shows the generic effect of a very small transverse 
perturbation on the otherwise regular longitudinal oscillation. As can be seen 
from the figures(3a) and (3b) that upto time = 600w p the velocity profiles 
remain same, but after that there is a catastrophic change in the velocity 
profile of V x and simultaneously the coherent oscillation in V z breaks up into 
a chaotic one. A fast Fourier transform of these oscillations had shown that 
the most dominant frequency for both V x and V z components are the same, 
suggesting energy equilibration with the onset of chaotic motion. The most 
intersting investigation here would be to see if this energy equilibration is 
also accompanied by color equilibration. There have been some effort in that 
direction by Gyulassy p7fl , but we will talk about it in the appropriate section. 
Having described the role of soft dynamics in RHIC, we'll now briefly introduce 
the corresponding picture in the still higher energy scale. That is the parton 
cascade model. As one moves to higher and higher energy, it's believed that the 
hard processes will dominate and produced partons will screen the long-range 
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field. Therefore to get the picture of partonic evolution in parton Cascade 
Model, one needs to modify equation (|l3|), by putting the external field to be 
equal to zero. Its worth noting that as a result of this, the two coupled partial 
differential equations decouples from each other and one is left with:, 

p»d ll W{x, P ,t)=C dUs ion (33) 

This collision term describes, all the collision-al effects due to 2 — > 2 and 
2 — > 1 and the reverse processes, with appropriate modification of distribution 
functions due to pauli blocking and bose enhancement, along with virtuality 
factors to take into account the space like and time like cascading. The patho- 
logical low pt divergence is cured by self consistently introducing a screening 
mass for the partons.This selfconsistently modifi ed p arton cascade model is 
known as the self screened parton cascade model pq j ,f or the sake of brevity 
we won't discuss it here, the details can be found m p3]. 



4.1 Discussion and outlook 

Having introduced the basic mechanisms and models of parton production 
in relativistic heavy ion collision, we will go over to the discussion of the 
prediction's of these models -valid at different energy scales. As before we 
start from the low energy non-perturbative soft process dominated regime. In 
the Boltzman Vlasov description with a Schwinger like source term, it appears 
that the initial temperature at the onset of equilibrium is dependent on the 
initial field strength that goes into the description. According to refrence Q for 

Eq = 20GeV/ fm 3 the initial temperature at the onset of thermalisation was 
found to be varying between 225 to 250 MeV. On the other hand according 
to reference ( [19 fig. a constituent quark of mass .5 GeV the initial peak 
temperature comes out to be around 4.5 Gev, though for current quark mass, 
the temperature comes out much less. According to reference ( [|30| ) the initial 
temperature in RHIC ( LHC) for gluons are estimated to be 500 MeV ( 660 
MeV) the same for quarks is 200 MeV ( 260 MeV) and the thermalisation time 
for Gluons varies from 0.3 to 0.25 fm/c in RHIC and LHC energy scales. For 
gluons and light quarks the predictions from [Q or |HJ see to be close to each 
other and they suggest that the system undergoes a multistage thermalisation. 
The quarks are believed to be less equilibrated than gluons [ 32 1. Although 
for charm or bottom quarks the thermalisation time ( 2 fm/c. The Initial 



Conditions as obtained from Self Screened Parton Model [33] for Gold on Gold 
at BNL RHIC and CERN LHC energies-prior to hydrodynamic expansion is 
given in Table I. 

Table I. 



Self Screened Parton Cascade 


Ti (fm/c) 


Ti 


e t {GeV/fm A ) 


RHIC 


0.25 


0.668 


61.4 


LHC 


0.25 


1.02 


425 
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Coming back to the question of color equilibration, though according to J27||, 
color equilibrates takes place before the system reaches thermal equilibra- 
tion, but it would be nice to demonstrate this result using the color kinetic 
equations, retaining their full nonlinearity. It is generally felt that such an 
investigation would shed light on the question, whether the onset of chaotic 
oscillation has really got some thing to do with the equilibration procedure 
or not. And lastly it's worth mentioning that all the previous studies reveal 
that in RHIC the system would not come chemical equilibrium though at LHC 
energy scale it might get very close to it. 
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